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Abstract 

We outline a general method for obtaining exact solutions of Schròdinger 
equations with a position dependent effective mass and compare the re- 
sults with those obtained within the frame of súper symmetric quantum 
theory. We observe that the distinct effective mass Hamiltonians proposed 
in the literature in fact describe exactly equivalent systems having iden- 
tical spectra and wave functions as far as exact solvability is concerned. 
This observation clarifies the Hamiltonian dependence of the band-offset 
ratio for quantum wells. 

PACS numbers: 03.65.Ca, 03.65.Ge, 03.65.-w 

1 Introduction 



An interesting question arises when one tries to formulate the correct Hamilto- 
nian for a particle with spatialiy dependent mass in an arbitrary potential well. 
This problem often arises in the calculation of confined energy states for carriers 
in semiconductor quantum well structures under the envelope-function and the 
effective-mass approximations where the effective mass of a carrier is spatialiy 
dependent on the graded composition of the semiconductor alloys used in the 
barrier and the well region of the nanostructures. Since the mass and the mo- 
mentum operators no longer commute, the correct ordering of these operators 
within the kinetic energy operator cannot be trivially assigned. This problem of 
ordering ambiguity is a long standing one in quantum mechanics, see for instance 
the excellent critical review by Shewell [1]. The study of quantum systems with 
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position dependent effective masses has been the subject of much activity in re- 
cent years. Apart from being an interesting tòpic itself such equations have found 
wide applications in the search of not only electrònic properties of semiconductors 
[2] but also quantum dots [3], quantum liquids [4], 3 He clusters [5], and metal 
clusters [6]. 

Although the treatment of effective-mass Schròdinger equations with non- 
constant mass is difhcult, some exactly solvable models for such systems have 
been recently introduced [7-11]. However, from the physics point of view, none of 
these works seem completely convincing. For example, Dekar and his co-workers 
[7] solved exactly the one-dimensional Schròdinger equation, derived from the 
general form of the effective-mass Hamiltonian without any restriction in the 
ordering parameters, for a system with smooth potential and mass step. The 
effective mass Hamiltonian and the connection rules for a system with abrupt 
heterojunction were deduced from the study of the limiting case when the smooth 
step potential and mass tend to an abrupt potential and mass step, which sug- 
gest that the appropriate form of the effective-mass Hamiltonian for an abrupt 
heterojunction is the form proposed by [12]. However, this modelling of the po- 
sition dependent potential and mass by a discontinuous profile is not physically 
acceptable, because in the real world neither the potential nor the effective mass 
can change abruptly across the heterojunction. In a different perspective, Du- 
tra and Almeida [8] discussed the relationship between exact solvability of the 
effective-mass Schròdinger equations and the ordering ambiguity. They intro- 
duced a simple technique to transform mass dependent Schròdinger equations 
to the conventional Schròdinger equation with a constant mass to obtain in a 
closed form of the full spectrum and wavefunctions of the original system. As an 
illustrative example they considered a physically plausible case (an exponentially 
changing mass and potential) and obtained the corresponding exact solutions by 
mapping the system onto the well known harmònic oscillator with centripetal 
barrier. It was nicely clarified that in this case the energy levels are redefined in 
such a way that the ordering ambiguity disappear. However their model has a 
drawback, which works only for a restricted choice of mass and potential functions 
to be able to serve exact solutions. Although the works in Refs. [9, 10] within 
the frame of supersymmetric quantum theory, together with a more recent work 
[11] used a Lie àlgebra, provide a powerful treatment for the systems of interest, 
all the formulations introduced in these models are only vàlid for the BenDaniel 
and Duke Hamiltonian [12] as they followed the work of Levy-Leblond [13]. 

Notwithstanding, taking into account the spatial variation of semiconduc- 
tor type, some effective Hamiltonians have been proposed with a coordinate- 
dependent mass for the carrier [12,14-18], and many authors have been trying 
to determine the correct Hamiltonian phenomenologically, for a review see [17]. 
Although, there has been growing consensus in favor of the BenDaniel-Duke 
Hamiltonian, see the related references in [17], the form of the effective-mass 
Hamiltonian has been still a controversial subject due to the location dependence 
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of the effective mass and the question of the exact form of the kinetic energy 
operator is still an open problem. In the light of all these, the present work ad- 
dresses a more general and compact framework introducing two different models 
to solve exactly for a physical system involving a position dependent mass with 
the consideration of all physically acceptable Hamiltonians proposed in the lit- 
erature. We will make clear that though the effective Hamiltonians appear in a 
different form, they describe in fact exactly equivalent systems having identical 
spectra and wavefunctions as far as exact solvability is concerned. Although the 
literature covered similar problems, to our knowledge an investigation such as 
the one we have discussed in this paper was missing. 

The paper is organized as follows. Using the spirit of the works [19] we intro- 
duce, in the next section, a generalized exact treatment procedure for effective- 
mass Hamiltonians and use the works in [9, 10] to extend their scenario and 
involve all possible Hamiltonians. Section 3 involves applications of the models 
developed for two different but physically meaningful mass considerations. Dis- 
cussion and analysis of the results obtained are given in section 4. Finally, some 
conlusions are drawn in section 5. 



2 Theoretical considerations 

The single-band effective-mass approach to a quantum-well problem requires that 
the envelope function ip satisfy the effective-mass equation 

Htp = Eip (1) 

where E is the energy eigenvalue and H is the Hamiltonian operator consisting 
of kinetic energy operator ( T ) and the potential energy operator (V"), 

H = f(z) + V(z) (2) 

Due to the compositional variation in a quantum well as a function of location, 
the kinetic energy and the potential energy are expressed as position-dependent 
operators in Eq. (2). The kinetic energy operator can be considered to be 
composed of four elements: 1 / \Jm(z), 1 / yjm(z), p, and p, where z and p are 



position and momentum operators, respectively. Because 1/ ym(z) and p are not 
commutable, there are different possible permutations to represent the kinetic 
energy operator [12,14-18]. All of these single-band effective-mass Hamiltonians 
are special cases of a general form of the Hamiltonian introduced by von Roos 
[20], 



H, 



1 



vR 



m a (z)pm í3 (z)pm' y (z) + m y (z)pm l3 (z)pm a (z) + V(z) , (3) 



where a + (3 + j = — 1. 
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By the correspondence in wave mechanics p — > —ihd/dz and z — > the 
effective mass equation, Eq. (1), togeher with any possible Hamiltonian proposed 
in the literature as the special cases of Eq. (3) can be written in a differential 
form, 

1 d-ip(z) 
m(z) dz 

where V e ^(z) is termed the effective potential energy whose algebraic form de- 
pends on the Hamiltonian employed (see Table 1), 

V eff (z) = V(z) + U aj (z) 



r^d_ 

' 2 dz 



+ V eff (z)4>(z) = EïP(z) 



(4) 



= V(z) 



fi 2 



[a + 7) m(z)m"(z) — 2 (cry + a + 7) m! (z) 



(5) 



4m 3 (z) 

in which the first and second derivatives of m(z) with respect to z are denoted by 
m! and m", respectively. The effective potential is the sum of the real potential 
profile V(z) and the modification U ay (z) emerged from the location dependence 
of the effective mass. A different Hamiltonian leads to a different modification 
term. 

In the following section we introduce two different theoretical models to solve 
Eq. (4) exactly, which lead us to define analytically solvable potentials for the 
systems undertaken. 



2.1 Coordinate transformation method 

In a recent work, Roy and Roy [19] have presented a coordinate transformation 
method to solve analytically an effective-mass Hamiltonian corresponding to the 
one proposed by BenDaniel-Duke [12]. We extend and generalize their model 
by tackling the problem with a fundamental point of view, ie., without using 
a particular form of the effective potential in (5) and arrive at a conceptually 
consistent result for exactly solvable systems with location dependent masses. 
Our generalization is based on the work of De et all [19] in which they gave 
explicit point canonical transformations and interrelated many shape-invariant 
exactly solvable potentials within the frame of the usual Schròdinger equation 
with a constant mass. We show that their simple approach, which consists of 
mapping through canonical transformations of coordinates, can also be used to 
interrelate two distinct systems having constant and position-dependent masses. 

As the aim is to solve Eq. (4) in a closed form to obtain the corresponding full 
energy spectrum and wavefunctions exactly for a solvable potential, one needs to 
transform Eq. (4) into a workable frame. Proceeding with a transformation of 
both the independent and dependent variables of the form 

z = f(z) , ip(z) = v(z)ip(z) , m(z) = m[f{z)\ = rh(z) , (6) 
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we transform Eq. (4) to 



2rh ] dz 2 



v 



f" m'\ d^J 
f m I dz 



v 



ml 



ell 

m j 



v 
v 



i/, 



+ fi{v*ff[f(z)}-E}4> = 0, (7) 



in which the prime denotes differentiation with respect to z . One can now easily 
reduce Eq. (7) to the Schròdinger equation with constant mass, which requires 



u(z) = Cy/f'(z)m(z) , f(z) = 



m 



rh(z) ' 



(8) 



where ítiq stands for a constant effective mass, and C is a constant of integration. 
Eq. (7) in this case reads 



~ 12 



1 m" 



2m dz 2 \ 2m \IQ m? 4 rh 



+ V e "[f(z)]\$(z)=E$(z) . (9) 



It is important to note that Eqs. (4) and (9) have identical spectra. With the 
wavefunctions ip(z) being square integrable, ie., (ip(z) \ip(z)) = 1 , setting the 
integral constant C = 1/^/mo i n Eqs- (6) and (8) and having in mind that 
m(z) > we find that (vp(z) \ip(z)) = 1 as well, ie., the functions tp(z) are also 
square integrable. 

For further simplification, we recali the mathematical definition used in ar- 
riving at Eq. (8), f' 2 (z) [rh(z)/mo] = 1, and remind that it may be inverted and 
written as 



= -L [Jmjz)dz = f-\z) , 



(10) 



which defines (though in implicit form) the mapping function f(z) and conse- 
quently enables finding rh(z) , together with the effective potential V e ^ [/(£)]■ 
Bearing in mind that 



m 



drh(z) dz d 



dz dz dz 
Eq. (9) can also be reduced to the form, 
h 2 d 2 ^j(z) f h 2 



m(z)] = f'(z)m'(z) 



m dm{z) 



m(z) dz 



2m dz 2 [32m 3 (z) 



7m' 2 {z) - Am(z)m"(z)] + V eff {z) $(z) = E^(~z) 

) z=z 

(12) 

It is worth to note that the change of variable in (10) may not always be invertible 
or not easily be invertible. But this does not really pose a problem as far as 
solvability of Eq. (4) is concerned. This is because the new coordinate z as a 
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function of the old coordinate z is explicitly known from (10) and if we choose 
V eff such that 

V eff (z) = V E s(z) - V m (z) , (13) 

where Ves is an exactly solvable potential and V m is the mass- dependent part of 
the full potential, 

fc2 



obviously one then finds exactly the spectrum of the system in (4). The corre- 
sponding wavefunctions can be obtained using Eq. (6). 

The interest in exactly solvable problems in quantum physics has increased 
sharply in the last few years. This is concerned of course with the fact that the 
behavior description of some physical systems is usually very complicated, but in 
some cases such systems can be modelled by means of a quite simple Hamiltonian 
which leads to Standard problems of quantum mechanics, as the one presented in 
this section. As the spectral properties of the constant-mass Schròdinger equation 
of solvable potentials are well known in the literature, one can readily obtain a 
corresponding potential for the effective-mass Schròdinger equation with identical 
spectra. The problem of generating isospectral potentials in quantum mechanics 
has been considered for more than 50 years, but recently the research efforts on 
this tòpic have been considerably intensified. 

In the next section, we discuss an alternative treatment for the problem under 
consideration and illustrative examples will be given in section 3. 



V„ 



7m' (z) — 4m(z)m"(z) 



2.2 Supersymmetric approach 



In the light of the recent works [9, 10], we introduce a superpotential W(z) and 
the associated pair of operators A and A + defined by 



A*(z) = . = -- ^ + W(z)1>(z) , A+4(z) = -± 
y/2m(z) dz dz 



+ W{z)i){z) 
(15) 



Notice that, due to the location dependence of the mass, d/dz and h/y2m(z) do 
not commute anymore. Within the frame of supersymmetric quantum mechanics 
[21], the partner Hamiltonians read 



H 1 = A + A = 



h 2 



d 



2m(z) dz 2 \2m(z) J dz 
2hW'(z) 



W 2 (z) - 



nwjz) 



H 2 = AA + = H 1 + 



^2m(z) ^2m(z) \^2m{í 



(16) 
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We see that the two partner Hamiltonians describe partides with the same 
efective mass-spatial dependence but in different potentials. As the usual exactly 
solvable potentials, as well as many recently discovered ones, are given in the form 
of the first partner potential in (16), we will work through the applications in 
the next section with the consideration that the effective potential in (5) should 
satisfy the condition, 



V eff ( Z ) =Vy = W 2 (z) - 




The second partner potential is used for obtaining the corresponding energy spec- 
trum. Considering the shape invariance concept [21], the supersymmetric partner 
potentials in (16) obey the integrability condition 

V 2 (z, ai) = Vi(z, a 2 ) + R{a x ) , a 2 = /(ai) . (18) 

The above relation means that the partner potentials V\ and V 2 have the same 
form, but are characterized by different vàlues of the parameter sets. The energy 

n 

eigenvalues, E n — J2 R( a i)> an d eigenfunctions of shape invariant poteantials 

i=i 

can be obtained in algebraic fashion. For more details the reader is refered to 
[9, 10, 21]. 



3 Applications 

A square quantum well subject no external field is usually modeled by a step 
profile which discontinues at the heterojunctions. The commonly used interface 
conditions at the heterojunctions are: (i) the continuity of the envelope function 
and its first derivative, or (ii) the continuity of the envelope function and its first 
derivative divided by the effective mass. Within each flat region of the square 
quantum well, the mass is a constant, and the effective potentials, Eq. (5), are 
identical because the derivatives of the mass with respect to the position vanish. 
Thus envelope functions within each flat region are independent of the Hamil- 
tonian used for the analysis. If the same interface conditions at the discontinue 
points between two adjacent regions are essentially imposed for the distinct ef- 
fective Hamiltonians, the eigensolutions will be exactly identical. As a result, the 
effective-mass Hamiltonians in (5) are expected to produce the same transition 
energy, namely, the band-offset ratio will be independent of the Hamiltonian if 
the heterojunction is modeled by a step function with essentially imposed inter- 
face conditions. This point will be discussed in the next section within the frame 
of our results. 

However, the discontinuity of the square quantum-well model implies an in- 
finite external elèctric field at the heterojunctions, and this is not physically 
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possible. In reality, the potential changes over a few monolayers for a perfect 
microscopic interface. In this section, the square quantum well is modeled by 
more realistic smoothed profiles, which leads to different effective potential de- 
scriptions. This modified realistic potential profiles remove the discontinuity of 
the sharp square quantum well. Considering the realistic quantum well applica- 
tions in the works [17, 22], we test our models introduced in section 2 employing 
the two different but physically meaningful effective mass variations discussed in 
the following. 



3.1 For m{z) = moexp(Az) 

To demonstrate the simplicity of the models used we consider here a particular 
case, as a first illustrative example, for which one have exact solutions. That is 
a particle with exponentially increasing or decaying, |À| oc 1/L with L being the 
quantum-well width, in the presence of a potential with similar behaviour that 
will be defined for each Hamiltonian in this section. 

Starting with Eqs. (12-14), within the framework of the coordinate trans- 
formation technique and taking the harmònic oscillator potential as a solvable 
potential, Ves(z) = Bz 2 , and bearing in mind that z = (2/À) exp (Xz/2) from 
Eq. (10), we obtain a variety of exactly solvable effective potential descriptions 
appearing in the original effective mass equation, Eq. (4). These are, correspond- 
ing to distinct effective mass Hamiltonians presented by Eqs. (3) and (5), 

V* f J D (z)=V eM\z)--—exp(-\z) , a = 1 = , (3 = -1 , (19) 

Ó2m 

for the BenDaniel-Duke Hamiltonian [12], 

hh 2 \ 2 

VèíUz) = V ex P (\z) + ^—ex P (-\z) , (3 = 7 = , a = -1 , (20) 
for the Gora and Williams (or the Bastard ) Hamiltonian [14, 15], and 

V^J K (z)=V eKp(\z) + ^-exp(-Xz) , a = 7 = -\ , = 0, (21) 

32m 2 

for the Zhu and Kroemer Hamiltonian [16]. We also obtain exactly the same 
result as in (21) for the Li-Kuhn Hamiltonian [17] (/3 — 7 = — 1/2 , a — 0). 
In the above equations, V = AB/\ 2 = h 2 \ 2 /32m . From Eqs. (4) and (12), 
it is clear that though the appearance and behaviour of the potentials in (19)- 
(21) are dissimilar, they have identical spectra, E n — (n + l/2)hu, where uj = 
\J2B /m . This result supports the similar discussion presented in the section 
3 of a recently published paper [11]. Moreover, from Eqs. (6) and (8), these 
Hamiltonians have the same wavefunctions, tp n (z) = N n [exp(Àz)] 1 ^ 4 ^ n {z) where 
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^ 1/2 2 

ï>n{z) = {2^7 #n (v^ 5 ) ex P ( _2 fr) with ^ = VZBmo/h. As a result, the 
distinct effective Hamiltonians considered here are not only isospectral but also 
describe identically equivalent systems as far as exact solvability is concerned. 
This significant result is also confirmed below with use of the supersymmetric 
approach described in section 3 as an alternative treatment. 

For the supersymmetric considerations, in the light of Eq. (35) of [9] we 
propose a superpotential 



w(2) = ^ exp(A , /2) __^ exp( _ A2/2) , (22) 

where 5 = hu , which yields the supersymmetric partner potentials in (16) in the 
form 

Vl(z) = V&exp(Az) - ^-^-exp(-Az) - ^ , 

3h 2 X 2 5 
V 2 (z) = V exp(Az) - ^— exp(-Xz) + - , (23) 

which is the simplest case of the shape invariance integrability condition given by 
(18) due to the partner potentials in (23) differing only by a uniform energy shift 
by 5 . Note that the first partner has the same shape as in (19) corresponding 
to the BenDaniel-Duke effective potential. Considering the shifting term 5/2 , 
together with Eq. (18), one can easily find the corresponding energy spectrum, 
E n = (n+l/2)hüü, which overlaps with the one found through the transformation 
technique. For the other exactly solvable shape invariant effective potentials we 
use a simple expression, 

V e JJ(z) = V e JJ D (z) - U ai {z) = Vï(z) - U ai {z) , (24) 

in which U ai {z) is the modification term in Eq. (5), 



U^ W (z) = exp(-A*) , Uf~ K {z) = U L a ~ K = exp(-Az) . (25) 



The substitution of Vi(z) in Eq. (23) into Eq. (24) leads us to arrive at 
Eqs. (20) and (21), which confirms the reliability of the coordinate transforma- 
tion technique developed in the previous section. Once more it is clear that the 
potentials in (19-21) will have the same energy spectra and wavefunctions as a 
unique superpotential is used, (Eq. 22), for the generation of these analytically 
solvable potentials within the frame of supersymmetric quantum theory. For the 
relation between the superpotential and wavefunction in case of the supersym- 
metric considerations the reader is referred to [9, 10]. 
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3.2 For m{z) = m (g^ 

This physically convenient mass variation, in which a is a positive constant and 
q(= Xz) involves the variable with a positive width parameter À oc 1/L, consid- 
ered here to convince the reader for that the models introduced in this paper 
works for all smoothly varying masses, unlike the recent models [7, 8]. There is 
another reason for the consideration of the work presented in this section. Among 
various potential shapes (quantum well profiles), there is one which has attracted 
some research attention recently: this is the case of Pòschl-Teller potential with 
a constant electron mass [23, 24] . Unfortunately, this idealized potential is not 
realizable in the common, graded ternary alloy based quantum wells, because of 
the effective mass therein necessarily varies, together with the potential. In this 
respect, the result shown in this section and its related discussion given in the 
next section would be helpful in designing realistic ternary alloy based structures 
with properties equivalent to those of idealized Pòschl-Teller potential. 

From Eq. (10), the relation between the transformed coordinate and the old 
one is z — z + (1/A)(a — l)ían _1 ç, and considering the Pòschl-Teller potential 
as one possible choice from analytically solvable potential family and working 
within the framework of the coordinate transformation method, one arrives at a 
class of exactly solvable potentials for the use in Eq. (4) belonging to different 
effective mass Hamiltonians, 

W/M-V m i ("-l)[3g 4 + g 2 (4-2a)-a]ft 2 A 2 
V eff (A v m (g ~ 1) [3g 4 + 2g 2 - a] ft 2 À 2 



Vl-M - V PT (z) + 2mo(g2 - fl)4 , 

V eff (z) _ v m (a-l)[3g 4 + g 2 (6-4a)-a]^ 2 A 2 

V BDD (z) - V PT (z) 2m (g2 + a) 4 ' ^ 

where V PT (z) — -A (A + ^) sech 2 {x[z+ (a ~ 1)f " n ' 1 q ] }. Note that for a -> 
1 ,m(z) — > m , all the above effective potentails reduce to the conventional 
Pòschl-Teller potential, V eff (z) — >• V PT (z), due to z -> z. From Eqs. (4), (8) 
and (12), the bound state energy spectra and wavefunctions corresponding to 
the potentials in (26) are E n = -(A- -^g=) and ip n (z) = N ny J*±£ rj){z) with 

ip{z) being the well described wavefunction [25] for the solution of the usual 
Schròdinger equation with a constant-mass Pòschl-Teller potential. 
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All these results are fully confirmed with use of the supersymmetric expres- 
sions presented in section 2 by adopting an ansatz for the superpotential, 



W(z) = A tanh \ \ 



z + 



l)tanh 1 q 



4 Discussion 

We discuss here the physics behind the results obtained. First, the applications 
given in the previous sections make clear the band-offset ratio dependence on the 
effective mass Hamiltonians, which is significant for quantum well applications. 
The conduction-band-offset ratio, which is the ratio of the conduction-band off- 
set to the total band gap of the heterojunction, has been investigated in various 
quantum wells because of its fundamental importance and application. The ratio 
has been measured by spectroscopic and electrical methods. From Duggan's and 
Kroemer's review articles [26] about the experimental and theoretical works, it 
can be seen that spectroscopic techniques are preferred over electrical ones in 
exploring the band-offset ratio via a quantum well, and that researchers using 
spectroscopy usually try to match their data with the theoretical results to de- 
termine the band offset. To demonstrate simply the band-offset ratio variation 
due to the choice of the Hamiltonian, we will focus here on transition energies 
between levels in conduction and valence bands. Following the works [17, 27] and 
considering only single-band effective-mass equations for the electron and the 
hole, Eq. (4), one finds the transition energy in the form of Et = E e + Eh + Eq 
where E e , E h are the eigenenergies in the conduction and valence bands, respec- 
tively, and E G is the band-gap energy. From the previous sections, it is obvious 
that the effective Hamiltonians undertaken will yield the same transition ener- 
gies between identical transition levels due to the identical vàlues for E e and 
Eh- Hence, the band-offset ratio for the BenDaniel-Duke Hamiltonian and for 
the others can be found by solving E T thers (Q others ) = E^ DD (Q BDD ), with Q be- 
ing the band-offset ratio. Recali the relation between the conduction band-offset 
ratio and the conduction-band-offset energy, V c b = QAE g , where AE g denotes 
the band-gap difference between binary and ternary materials. Having in mind 
that the band-offset ratio of a quantum well determines the barrier height of the 
conduction band and valence bands, V c b corresponds to the depth of the effec- 
tive potentials discussed through the article. This leads to the realization of the 
band-offset ratio dependence on the effective-mass Hamiltonian due to the un- 
derlying differences between the strengths of the potentials obtained. As a result 
of this, in the interpretation of a given spectrum, the Hamiltonian employed in 
the analysis cannot be regarded independently of the band-offset ratio utilized, 
unlike the case of a simple square quantum well consideration. Therefore, an 
attempt to determine the band-offset ratio from experimental data by matching 
calculated transition energies with spectral peaks would involve large inaccuracy. 
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There is another interesting point behind the present results. A systematic 
procedure based on the inverse spectral theory and supersymmetric transfor- 
mations has been recently proposed [24] for optimized design of semiconductor 
quantum well structures via tailoring the quantum well potential, which enables 
shifting bound states in a quantum well and makes the search for the best desired 
energy spectrum and potential shape. By varying the free parameters appearing 
in the procedure one can then design a convenient optimized structure. However, 
in these notable works the idealized constant-mass Pòschl-Teller potential, which 
allows one to set analytically the spacing between states, considered since a direct 
implementation of a more realistic position-dependent effective mass related to 
the position-dependent potential in their theoretical considerations is not trivial. 
The optimization of continuously graded structures thus require more sophisti- 
cated techniques. In this respect, we believe that the applications given in the 
previos section, in particular the one involved the Pòschl-Teller potential, give 
a lot of material for experimenting in the optimized quantum well làser design. 
For instance, though we have explicitly shown that the energy spectra of the 
realistic ternary alloy based structures with a carrier having a spatially varying 
mass are equivalent to those of the constant-mass Pòschl-Teller Potential, the 
maximization of the gain may be accomplished via changing the quantum well 
profile which is in turn changes the wave functions. Hence the consideration 
of our results within supersymmetric transformations, as an alternative to the 
recent procedure proposed, in particular the significant difference between the 
wavefunctions by v(z) shown in (6), relating the solutions of the Schròdinger 
equations with constant and location dependent masses, and use of the more 
realistic effective potentials in (26) in tailoring process instead of the Standard 
Pòschl-Teller potential, should reproduce more reliable results in order to help 
for the best design of such structures. Along this line the work is in progress. 

Moreover, Plastino and his co-workers [28] recently studied some simple one- 
dimensional quantum mechanical systems characterized by a piecewice flat po- 
tential and mass to illustrate the influence of a non-constant mass on the density 
of the bound state energy levels. With the consideration of a finite potential 
well they showed that the number of bound states is less than those of the con- 
stant mass situation when the effective mass inside the well is lower than that of 
outside (m ), and the opposite behaviour occurs when the effective mass inside 
the well is larger than the mass outside. However, our applications in section 
3 do not confirm their work. This contradiction may raise a further discussion 
on the reliability of the present results, which can be clarified as follows. In the 
two different variable mass definitions, m(z), for the quantum well used in the 
previous section, m(z) > itlq for the case À > and a > 1 while m(z) < itiq 
in case À < and a < 1. The consideration of m(z) > m case leads to sin- 
gle potential wells, sharper than that of the Standard potential corresponding to 
constant-mass potentials, whereas m(z) < m Q case give rise to bistable double 
well potentials, like the related illustrations in [9]. In spite of the different aspects 
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exhibited by the effective potentials defined in section 3, we have clearly shown 
that they share the same energy spectra regardless of À and a vàlues, together 
with the constant-mass potentials (À = 0, a = 1). Consequently, one expects that 
the density of bound states for the systems we interest should not depend on the 
variation of a carrier mass. 

5 Concluding remarks 

In this article, we have discussed the problem of solvability and ordering am- 
biguity in quantum mechanics as the form of the effective-mass Hamiltonian 
has been a controversial subject due to the location dependence of the effective 
mass. It was shown through particular examples that the exact solvability de- 
pends not only on the form of the potential, but also on the spatial dependence 
on the mass. Within the framework of the two different theoretical treatments, 
the effective-mass Schròdinger equation has been transformed to a constant-mass 
Schròdinger equation and we have clarified that the Schròdinger equations with 
different masses and potentials can be exactly isospectral. We have also shown 
that though the potential energy of the BenDaniel-Duke Hamiltonian differs from 
the effective potentials of other Hamiltonians proposed in the literature by a term 
caused by the mass dependence on location, the exact analytical solutions to the 
effective-mass equations do not change with the Hamiltonian. As far as we know, 
this feature was not perceived until now. The discussion given behind the results 
obtained may be of interest, e.g., in the design and optimization of semiconductor 
quantum wells. In addition to their practical applications, we believe that the 
study of quantum mechanical systems with a position-dependent mass within the 
framework of the present models will raise many interesting conceptual problems 
of fundamental nature. In particular, the methods used in this article may be 
extended to find applications in also the study of quasi- and conditionally-exactly 
solvable systems with non-constant masses. 



13 



Table 1: Single-band effective mass Hamiltonians, Eq. (3). 



Hamiltonian 


a 
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Ref. [12] 





-1 





Ref. [14, 15] 


-1 








Ref. [16] 


-1/2 





-1/2 


Ref. [17] 





-1/2 


-1/2 



14 



References 

[1] J. R. Shewell, Am. J. Phys. 27, 16 (1959). 

[2] G. Bastard, Wave Mechanics Applied to Semiconductor Heterostructures 
(Les Editions de Physique, Les Ulis, France, 1992). 

[3] L. I. Serra and E. Lipparini, Europhys. Lett. 40, 667 (1997). 

[4] F. Arias de Saavedra, J. Boronat, A. Polls and A. Fabrocini, Phys. Rev. B 
50, 4248 (1994). 

[5] M. Barranco, M. Pi, S. M. Gatica, E. S. Hernandez and J. Navarro, Phys. 
Rev. B 56, 8997 (1997). 

[6] A. Puente, Li. Serrea, and M. Casas, Z. Phys. D 31, 283 (1994). 

[7] L. Dekar, L. Chetouani and T. H. Hammann, Phys. Rev. A 59, 107 (1999); 
J. Math. Phys. 39, 2551 (1998) 

[8] A. De Souza Dutra, C. A. S. Almeida, Phys. Lett. A 275, 25 (2000). 

[9] A. R. Plastino, A. Rigo, M. Casas, F. Gracias and A. Plastino, Phys. Rev. 
A 60, 4398 (1999). 

[10] V. Milanovic and Z. Ikovic, J. Phys. A : Math. Gen. 32, 7001 (1999). 

[11] B. Roy and P. Roy, J. Phys. A : Math. Gen. 35, 3961 (2002). 

[12] D. J. BenDaniel, C. B. Duke, Phys. Rev. B 152, 683 (1966). 

[13] J. M. Levy-Leblond, Phys. Rev. A 52, 1845 (1995). 

[14] T. Gora, F. Williams, Phys. Rev. 177, 1179 (1969). 

[15] G. Bastard, Phys. Rev. B 24, 5693 (1981). 

[16] Q. G. Zhu, H. Kroemer, Phys. Rev. B 27, 3519 (1983). 

[17] T. Li, K. J. Kuhn, Phys. Rev. B 47, 12760 (1993). 

[18] F. S. A. Cavalcante, R. N. Costa Filho, J. Ribeiro Filho, C. A. S. De Almeida 
and V. N. Freire, Phys. Rev. B 55, 1326 (1997). 

[19] R. De, R. Dutt and U. Sukhatme, J. Phys. A : Math. Gen. 25, L843 (1992); 
B. Roy and P. Roy, quant-ph/0106028. 

15 



[20] O. Von Roos, Phys. Rev. B 27,7547 (1983). 

[21] F. Cooper, A. Khare and U. Sukhatme, Phys. Rep. 251, 267 (1995). 

[22] S. Datta, Quantum Phenomena, Modular Series on Sòlid State Devices, Vol. 
8 (Addison-Wesley, Reading, Ma, 1989), Chaps. 1 and 6; S. Adachi, J. Appl. 
Phys. 58, Rl (1985); P. Lawaetz, Phys. Rev. B 4, 3460 (1971). 

[23] S. B. Haley, Am. J. Phys. 65, 237 (1997); B. Y. Tong, Sòlid State Commun, 
104, 679 (1997); B. Y. Tong, N. Kiruischeva, Phys. Lett. A 29, 49 (1997). 

[24] G. Todorovic, V. Milanovic, Z. Ikonic and D. Indjin, Phys. Lett. A 279, 268 
(2001); S. Tomic, V. Milanovic, Z. Ikonic, IEEE J. Quant. Elect., 37, 1337 
(2001); S. Tomic, V. Milanovic, M. Tadic and Z. Ikonic, Sòlid State Comm., 
113, 221 (2000); S. Tomic, V. Milanovic, Z. Ikonic, Òptics Comm., 143, 214 
(1997). 

[25] S. Flügge, Practical Quantum Mechanics, Springer-Verlag, Berlin, 1971. 

[26] G. Duggan, J. Vac. Sci. Tech. B 3, 1224 (1985); H. Kroemer, Surf. Sci. 174, 
299 (1986). 

[27] A. Chomette, B. Deveaud, M. Baudet, P. Auvray and A. Regreny, J. Appl. 
Phys. 59, 3835 (1986). 

[28] A. R. Plastino, A. Puente, M. Casas, F. Garcias and A. Plastino, Revista 
Mexicana De Fisica. 46, 78 (2000). 



16 



